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We estimate the quark-mass dependence of the topological susceptibility with dynamical overlap and clover 
fermions. Unquenching effects on the susceptibility turn out to be well approximated by a reweighting of a 
quenched ensemble with a low-eigenmode truncation of the fermionic determinant. We find that it is most likely 
due to the explicit chiral symmetry breaking of the fermion action that present day dynamical simulations do not 
show the expected suppression of the topological susceptibility. 



One of the most profound effects the inclusion 
of dynamical fermions is expected to have on the 
QCD vacuum is the suppression of fluctuations of 
the topological charge. Due to the index theorem 
in the continuum, the Dirac operator has \Q\ zero 
modes in the charge Q sector. In the presence of 
Nf degenerate flavours of light dynamical quarks, 
in the charge Q sector the fermionic determinant 
is proportional to m^ Nf (to lowest order in m). 
Light dynamical fermions are therefore expected 
to suppress higher topological sectors and thus 
also the topological susceptibility, \ = (Q 2 )/V- 

On a more quantitative level, for small enough 
quark masses chiral perturbation theory predicts 
that 
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where S is the chiral condensate and m n and 
/„• are the pion mass and decay constant [Q. 
This holds provided the volume is large enough 
that chiral symmetry is effectively broken, i.e. 
x = T,mV 1. 

At the other extreme, i.e. for heavy quarks, 
one expects the susceptibility to approach a con- 
stant, namely the quenched value x quonch cd = 
(203 ± 5MeV) 4 g). The interpolation of the sus- 
ceptibility from light to heavy quarks has recently 
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been discussed by S. Durr ||. He noticed that 
besides the light dynamical quarks the finite vol- 
ume of the box also suppresses fluctuations of the 
topological charge, since only a finite number of 
instantons can be accomodated in a finite box. 
Assuming that the two suppression mechanisms 
work independently, Durr derived and essentially 
unique formula for the susceptibility. In terms of 
the pion mass it reads 
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having two free parameters, and X auMlolled . 
Durr also pointed out that — most probably due 
to scaling/chirality violations — all the presently 
available dynamical susceptibility data lie signif- 
icantly above this curve. 

We expect that the suppression of the suscep- 
tibility depends crucially on the index theorem 
which in general does not hold on the lattice 
The reason for this is that Dirac operators that 
explicitly break chiral symmetry (e.g. the Wilson 
operator) do not have exact zero modes and be- 
sides, the very definition of the topological charge 
is ambiguous on the lattice. It is thus not clear 
how exactly light dynamical fermions suppress 
higher |Q| sectors on the lattice. The aim of the 
present work is to check whether improving the 
chirality of the fermion action can have a signif- 
icant effect on the topological susceptibility. In 
what follows, the topological charge is always de- 
fined with my overlap Dirac operator by counting 
zero modes. Note that this implies an exact "in- 
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dex theorem" which in fact is a tautology with 
this charge definition. 

Unfortunately chiral fermion actions are still 
tremendously expensive and a full-fledged dy- 
namical simulation is out of question. It turns out 
however that in the present case that might not 
be needed. A reasonable estimate of the quark- 
mass dependence of the susceptibility can be ob- 
tained by reweighting a quenched ensemble with 
a low-mode truncation of the fcrmionic determi- 
nant. Computing the lowest N eigenvalues of the 
Dirac operator D, the determinant can be written 
as 

N 

det(£> + m) « JJ(A fc + m) x (3) 

fe=i 

which in turn provides the full quark-mass and 
Nf dependence of the topological susceptibility. 
Notice that while the truncation (||) with N small, 
does not provide a good approximation to the 
full determinant, the topological charge is corre- 
lated only with the low-end of the Dirac spec- 
trum. Therefore, the difference in the effective 
action between the various charge sectors can be 
reliably estimated even with N reasonably small. 

Reweighting is generally a bad procedure since 
control over the errors is lost exponentially fast 
with increasing volume. In certain cases, however 
it can happen that there is still a large enough 
window in the volume where interesting results 
can be obtained (see e.g. || for a recent exam- 
ple). In the following I give some reasons why 
reweighting could work for the topological sus- 
ceptibility. 

First of all, there are two sources of the fluc- 
tuations in the fermionic determinant. Eigenval- 
ues A/c with k > N, which is by far the largest 
contribution since with our choice of N, most of 
the eigenvalues belong to this sector. This part 
does not appear in our approximation. For light 
quarks an important contribution to the fluctua- 
tions comes from the zero modes but this is not 
harmful since this is exactly the quantity we want 
to measure: the difference in the effective action 
between the various charge sectors. Secondly, as 
the sea quark-mass changes, in principle, all the 
other parameters of the theory (lattice spacing, 
etc.) change. But this is not the case if only 



the low eigenvalues are used for the reweighting. 
Clearly e.g. the heavy quark potential is expected 
to be the same in all topological sectors. In spite 
of all this, the results that I present here have to 
be considered as preliminary. The error bars are 
estimated with a jack-knife analysis and a more 
careful error analysis will appear in a future pub- 
lication. This is certainly needed to be able to 
judge how far one can trust the results in the 
small quark mass limit. 

The simulations were done with the Wilson pla- 
quette action at (3 = 5.85 corresponding to a lat- 
tice spacing of about 0.12 — 0.13fm. The volume 
of the box was set to 2.5fm 4 which turned out to 
be a good compromise; this is already a physically 
reasonable volume (if somewhat small) and at the 
same time the reweighting still works reasonably 
well. The overlap action that I use is constructed 
from 10 times APE smeared fat links || and the 
c sw = 1-0 clover action. The results are based 
on 300 configurations and the full run took about 
10-20 PC-months. 
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Figure 1. The spectral density of the overlap 
Dirac operator in different Q sectors. 

The first question that arises is where exactly 
to truncate the determinant. To find this out, 
in Fig. [j] I show the spectrum of the overlap 
Dirac operator in the different topological charge 
sectors. There does not appear to be a differ- 
ence between the various Q sectors above roughly 
A > 0.20. Therefore I truncate the determinant 
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around A = 0.2 — 0.4 which at the given volume 
translates into taking the lowest N — 10 — 30 
eigenvalues. 
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Figure 2. The unqucnchcd topological suscepti- 
bility vs. pion mass with overlap fermions. 

In Fig. H I show the pion mass dependence 
of the topological susceptibility computed with 
truncations of N = 16, 24, 32. The steep line 
starting at the origin is the Leutwyler-Smilga pre- 
diction ([!]) and Durr's formula is also included us- 
ing U = 93MeV and x quenched = (200MeV) 4 . The 
horizontal shaded region is the quenched estimate 
for the susceptibility from the same data set. One 
expects the predictions based on broken chiral 
symmetry to be valid only if x 3> 1. At the given 
volume the value of the parameter x happens to 
roughly coincide with (m^ro) 2 , so chiral symme- 
try is effectively broken if (m^ro) 2 3> 1. There 
are some indications that even in that regime, 
there are considerable finite volume corrections 
to the unquenched susceptibility, much more so 
than in the quenched case. A more detailed study 
of this question will appear in a separate publica- 
tion. 

The same method can be used to compute the 
susceptibility with the clover action. For this I 
take the non-perturbatively determined clover co- 
efficient c sw = 1.91 0. The results are shown 
in Fig. H in a fashion similar to Fig. ||. Here I 
also plotted the UKQCD results || obtained with 
Nf — 2 0(a) improved Wilson fermions. 
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Figure 3. The unqucnchcd topological suscepti- 
bility vs. pion mass with clover fermions. 



On the same set of configurations the clover 
action produces a significantly higher susceptibil- 
ity than the overlap. This strongly indicates that 
chirality of the fermion action can be important 
for the topology. The error bars in the clover 
case are bigger than for the overlap, mainly due 
to the weaker correlation between the topological 
sectors and the fermionic determinant. This is 
even more pronounced in the case of unimproved 
Wilson fermions (not shown in the Fig.), where 
this makes it impossible to obtain sensible results 
using reweighting and the available statistics. 
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